The ground state of binary systems with a periodic modulation of the linear coupling 
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We consider a quasi-one-dimensional two-component systm, described by a pair of Nonlinear 
Schrodinger/Gross-Pitaevskii Equations (NLSEs/GPEs), which are coupled by the linear mixing, 
with local strength Q, and by the nonlinear incoherent interaction. We assume the self-repulsive 
nonlinearity in both components, and include effects of a harmonic trapping potential. The model 
may be realized in terms of periodically modulated slab waveguides in nonlinear optics, and in Bose- 
Einstein condensates too. Depending on the strengths of the linear and nonlinear couplings between 
the components, the ground states (GSs) in such binary systems may be symmetric or asymmetric. 
In this work, we introduce a periodic spatial modulation of the linear coupling, making Q an odd, 
or even function of the coordinate. The sign flips of Q(x) strongly modify the structure of the GS in 
the binary system, as the relative sign of its components tends to lock to the local sign of Q. Using 
a systematic numerical analysis, and an analytical approximation, we demonstrate that the GS of 
the trapped system contains one or several kinks (dark solitons) in one component, while the other 
component does not change its sign. Final results are presented in the form of maps showing the 
number of kinks in the GS as a function of the system's parameters, with the odd/even modulation 
function giving rise to the odd/even number of the kinks. The modulation of Q(x) also produces a 
strong effect on the transition between states with nearly equal and strongly unequal amplitudes of 
the two components. 

PACS numbers: 03.75.Lm; 42.65.Tg; 42.65.Wi 

I. INTRODUCTION 

Currently available experimental techniques make it possible to create two-component Bose-Einstein condensates 
(BEC). In the first experiments on such systems mixtures of different hyperfine atomic states were used, such as atoms 
of 87 Rb with different values of the total spin, F [1] . Such binary BEC systems offer straightforward possibilities for the 
realization of novel quantum phase transitions. In particular, while the superfluids formed by 87 Rb atoms in states 
\F = l,mp = —1) and \F = 2,mp = +1) are immiscible, the proximity of this system to the miscibility threshold 
makes it possible to create specific phase-separation patterns [3]. On the other hand, the pair of \F = l,mp = 1) 
and \F = 2,rriF = 2) states form a miscible supcrfluid, in which the separation between the species may be induced 
by an external field or by a specific flow pattern [3]. The immiscibility- miscibility border in binary systems may be 
shifted by means of the Feshbach-resonance technique, which changes the strength of the nonlinear interactions using 
external magnetic or optical fields [5] . 

Novel avenues in the physics of multicomponent BECs were opened with an realization of all-optical trapping and 
creation of spinor condensates, where the different components correspond to the different Zeeman levels of a given 
hyperfine manifold of spin F (for reviews cf. [6-8 ]). In this case , in principle, the number of (relevant) components 
can be controlled (from 2 to 2F + 1) using optical methods, and Feshbach resonances may be used to tune the desired 
collisional channels. 

Very accurate models for the description of BEC in quantum gases are provided by the Gross-Pitaevskii equation 
(GPE), or systems of coupled GPEs, in the case of the binary system [2]. With the help of these equations, ground 
states (GSs) of the binary systems [§] and excitations on top of them [TU] have been studied in detail. A well-known 
pattern created by the immiscibility in the binary BEC is the domain wall between spatial regions occupied by the 
different species [IT], which may also be accurately described by means of coupled GPEs |12j . 

An additional physically interesting ingredient of the binary-condensate settings is a possibility to induce the 
interconversion in the mixture of two hyperfine, or Zeeman states of the same atom by an electromagnetic spin-flipping 
field [13j . or a two-photon Raman transition (cf. [6]). In terms of the coupled GPEs, the interconversion is accounted 
for by linear-coupling terms. Diverse effects have been predicted in such models, including Josephson oscillations 
between the states [T3] , domain walls [THj , "co-breathing" oscillation modes in the mixture [TB] , nontopological vortices 
|17j , a linear-mixing- induced shift of the miscibility-immiscibility transition |18j , spontaneous symmetry breaking in 
two-component gap solitons [19], etc. 

Similar models, combining the nonlinear and linear interaction between two modes, occur in nonlinear optics. A well- 
known example is the system of coupled nonlinear Schrodinger equations (NLSEs) for two orthogonal polarizations of 
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light in a nonlinear optical fiber. In that case, the linear coupling represents the mixing between linear polarizations 
in a twisted fiber, or between circular polarizations in a fiber with an elliptically deformed core |20j . To develop 
a model in optics similar to that describing the BEC, one can consider a slab waveguide with a built-in array of 
fiber-like elements (strips or photonic nanowires buried into the slab, which gives rise to the linear mixing of TE and 
TM-polarized modes [5T]). A more versatile optical setting may be based on slab waveguides made as a photonic 
crystal [22]. Then, the propagation distance in the respective NLSE is a counterpart of the temporal variable in the 
GPE, while the transverse coordinate in the waveguide model plays the same role as the spatial coordinate in the 
GPE. 

An issue of obvious physical interest is the identification of the GS in linearly-coupled binary systems. A general 
limitation on all stationary states in such systems is that both components must have equal values of the chemi- 
cal potential [this condition is not imposed onto systems with the incoherent, i.e., XPM (cross-phase-modulation) 
interaction between the components]. Nevertheless, the GS in the symmetric linearly coupled system may feature 
a nontrivial asymmetric shape. Typically, the spontaneous symmetry breaking of the GS happens when the linear 
coupling competes with the nonlinear self-attraction acting in each component [23] . 

The relative sign of the two stationary fields is a significant characteristic of the GS in the linearly-coupled binary 
system, as this sign couples to that of linear-mixing constant (f2) through the corresponding term of the Hamiltonian 
density, hence it affects the choice of the state providing for the minimum of the Hamiltonian (see more details 
below). This circumstance is expected to become nontrivial if f2 is subject to a spatial modulation, which makes it 
a sign-changing function of the coordinate, Q = fi(x). Indeed, if the local structure of the GS adiabatically follows 
the change of f2(x), this means that the relative sign of the two components, being locked to the sign of fi(x), must 
flip along with it. Therefore, each zero crossing of Q(x) must give rise to a kink, alias dark soliton, in one of the two 
components, while the other one keeps a constant sign. This idea has been developed in Ref. [Mj, where the authors 
designed a method to generate dark solitons and vortices in multicomponent BECs by means of stimulated Raman 
adiabatic passage (STIRAP) from one internal state to another. While the initial BEC had a ground state wave 
function (without sign changes, or topological defects), the target BEC exhibited dark solitons or vortices depending 
on the spatial dependence and symmetry of the Raman coupling. 

The objective of this work is develop the ideas of Ref. [23] further and identify in more detailed and systematic 
way the structures of the GS in the linearly-coupled binary systems with f2(x) subject to the periodic change of its 
sign. This situation can be readily implemented in the above-mentioned optical model, by means of an appropriate 
superstructure created on top of the slab carrying the periodic waveguiding array or the photonic-crystal structure. In 
terms of the BEC, the periodic modulation of f2(x) may be realized too, in principle, if the spin-flipping electromagnetic 
field is patterned as a standing wave. Obviously, the latter realization cannot be achieved using a direct (microwave) 
transition, as the frequency of the transition between the hyperfine atomic states, ~ 10 GHz citeBBS, corresponds 
to wavelength ~ 30 mm of the electromagnetic field. This is too large in comparison with the size of experimental 
setups, and it could at best allow for a linear in space coupling, that could generate a single soliton, as pointed out 
in Ref. [24]. Note, however, that periodic modulation of f2(x) of desired period can be easily realized, if we use the 
two photon Raman transitions. 

It is also worth noticing that effect discussed here have a lot of similarity to the effect of disordered-induced-order 
[25j . which occurs in systems with continuous symmetry (such as complex phase, U(l) symmetry in BECs), when a 
symmetry breaking random, quasi-periodic, or even periodic field (coupled to the order parameter) is applied. Such 
random field will force the systems to order in a direction "orthogonal" to the symmetry breaking field. For example, 
a two-component BEC with real random (pseudo-random) Raman coupling between the components, will order in 
such a way that the relative phase between the components will be tt/2 [2B] . Similar, phase control effects and disorder 
induced orderings occur in Fermi superfluids [27] . or in quantum XY spin chains in ID [28] , 

This paper is organized as follows. The model, based on the system of coupled NLSEs/GPEs, is formulated in 
Section II, where some analytical results are reported too. In addition to the linear and XPM couplings, the equations 
also include a weak trapping potential, to make the extension of the GS effectively finite. Basic numerical findings, 
in the form of actual GS profiles and maps summarizing the most essential characteristic of the GS, viz., the number 
of kinks in it, are reported in Section III. The paper is concluded by Section IV. 
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II. THE MODEL AND ANALYTICAL RESULTS 



A. The coupled equations 



We consider a system of equations for wave functions of the two components, "01,2: which are cast in the notation 
corresponding to scaled GPEs, 



- l -V 2 + V (x,y) + \^ n \ 2 + g\^ n \ 2 
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with n = 1,2. The normalization is used to set the SPM (self-phase-modulation) coefficients equal to 1, while 
g is the free XPM coefficient. Equations (fTl) are written in the two-dimensional form, but assuming a strongly 
anisotropic trapping potential, V (x, y) = (1/2) (uj x x 2 + ^ y y 2 ) , where uj x — O.OlWj,, and ui x = ui = 0.05 is fixed below. 
Accordingly, the equations were solved numerically in a strongly anisotropic domain, —25 < x < +25, —0.5 < y < 
+0.5. The modulated linear-coupling constant was taken in two forms, odd and even ones: 



fl(x) — A {sin (ax) , cos (ax)} 



(2) 



(most results are reported below for the former shape of the modulation). The amplitude and wavenumber of the 
modulation, A and a, will be varied below as free parameters. As concerns the XPM coefficient, numerical results 
are reported below for g = —1, 0, 2/3, 1, and 2. In terms of the optics, the most relevant values are g — 2/3 and g = 2, 
which correspond, respectively, to the pairs of linearly and circularly polarized waves [2D]. The intermediate value 
g = 1 and "exotic" ones, g = and g = — 1 are, in principle, possible in photonic-crystal slabs. In the BEC system, 
both positive and negative values of g may be adjusted by means of the Feshbach-resonance method. 

Stationary solutions to Eqs. ([T| are sought for in the customary form, V+2 (x,y,t) = exp (— i/it) cf>i 2 (x, y). Real 
stationary wave functions (f>i,2 (recall chemical potential /i is identical for both components) were found by means of the 
standard numerical technique based on the solution of Eqs. ([I]) in imaginary time [29] . For analytical considerations, 
we adopt the stationary equations in the one-dimensional form, 



with fl(x) taken as per Eq. Q. 
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B. Analytical considerations 

To understand basic properties of the GS predicted by Eqs. ([3]), we start with the case of the uniform system 
(infinitely long or subject to periodic boundary conditions), i.e., with uj — 0, and a — 0, the linear-coupling coefficient 
corresponding to cos (ax) = 1 in Eq. ^ . The Hamiltonian density of the uniform states is 



H = 



+ <t>i)+g<t>l<i>l + A<t>i<t>. 



(4) 



A natural objective is to identify the GS as a state that minimizes the Hamiltonian density Q for a fixed value of 
the total atomic density, in terms of the BEC (or total-power density, in terms of optics) , Af = cf) 2 + 4> 2 . 

One can find two different uniform solutions to Eqs. ^ with lo = a = 0, symmetric and asymmetric. The former 

(5) 



one is 

0! = -sgn(A)0 2 , 4>\ = $ = N/2, n = (1 + g)Af/2 - \A\ /2, 
with the respective Hamiltonian density 



U 
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s ymm = ^(5 + 1) AA 2 --M\A\ 
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[a symmetric solution subject to the opposite sign locking, (f>i — +sgn(A)(/>2, exists too, but it gives rise to a higher 
Hamiltonian density]. The asymmetric solution is 



Af ± y H 2 - (g - 1) A 2 



(7) 
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(+ and — pertain to 4>\ and <f>2 ), with the relative sign of the wave functions locked so as to satisfy sgn((/>i</>2) = 
— sgn((g — 1) A). The corresponding values of the Hamiltonian density is 

1 I A 2 

Hasyxmu = -J\ - - - — ^ • (8) 

Obviously, asymmetric solution Q exists under the condition of \g — 1| > |A|/A/", and the comparison of expressions 
^ and ([8]) demonstrates that the GS, which must realize the minimum of the Hamiltonian density, corresponds 
to the asymmetric state at g > 1, and to the symmetric one at g < 1. This result exactly coincides with the 
elementary immiscibility (miscibility) condition in the uniform medium, g > 1 (g < 1) |30| . Thus, in the ideal uniform 
configuration, the point of the transition between the symmetric and asymmetric states, g = 1, is not shifted by the 
linear mixing. However, the difference from the linearly uncoupled system is that, in the immiscible phase (at g > 1), 
the binary system without the linear coupling must form domain walls , while the linearly coupled system may stay 
uniform, spontaneously concentrating a larger number of atoms (or larger optical power, in terms of optics) in one 
component. 

The state with a strongly broken symmetry, i.e., \4>2\ <C \<fii\, which corresponds to small values of \A\, can be found 
as an approximate nonuniform analytical solution to the full system of equations (|3), which include the trapping 
potential ( to > 0), and a^O. In the zero-order approximation (A — 0), one may take the solution as 02 = and 
4>i (x) in the form of the Thomas- Fermi (TF) ansatz. After that, the first-order (~ A) solution for 4>2(x), and the 
solution for <fii(x), which includes the second-order (~ A 2 ) correction, are found in the following form: at x 2 < 2p/uj 2 , 



, , , A v / p-u} 2 x 2 /2 J sin(ax) , 



2 (g - 1) (p - uj 2 x 2 /2) + a 2 /2 \ cos(aa; 



v {x) = y/p — oj 2 x 2 /2 s 1 T 



.4 2 



cos (2ctx) 



16 (a 2 + p - uj 2 x 2 /2) [{g -l)(p- lo 2 x 2 /2) + a 2 /2] 



(10) 



and (fix = <p2 = at x > 2p/oj . Here, \/ p — lo 2 x 2 /2 is the usual form of the TF approximation for 0i(x), with the 
effective chemical potential including a small correction, p = /i + A 2 / (4a 2 ), which eliminates a formal secular term 
at the second order of the perturbative expansion. The upper and lower signs in front of A 2 in Eq. ( 10 ) correspond, 
respectively, to the upper and lower rows in Eqs. ^ and Q. 



III. NUMERICAL RESULTS 



A. Ground-state profiles 

Simulations of Eqs. ([I]) in imaginary time t = —it were performed in the interval of < r < 500, with time step 
At = 10~ 3 . In fact, the convergence to stationary patterns was achieved by r ~ 100 — 200. The finally generated 
patterns do not depend on the choice of the seed configuration in the imaginary-time integration. 

A set of typical profiles of the GSs generated by the numerical method are displayed in Fig. [l] The profiles were 
obtained varying the amplitude of the odd modulation function in Eq. ([3]), A, and the XPM coefficient, g, for a 
particular fixed value of the modulation wavenumber, a = 0.4624 (other values a < 1 produce essentially the same 
picture, see also below). 

The approximate analytical solution given by Eqs. (J9J, ( [iO] ) adequately describes the general shape of the patterns, 
including the facts that the number of kinks is odd in the case of the odd modulation [and even, in the case of the 
even modulation, as shown by the numerical results generated for function cos (ax) in Eq. ([3])]. The dip in the profile 
of the stationary wave function <fii(x), which is observed, around x — 0, at g = — 1 for A = 0.01, and at g = —1, 0, 2/3 
for A — 0.1, is correctly explained by Eq. ( [To] ) , which includes the correction to the TF (zero-order) approximation 
with the curvature at the center opposite to that of the TF waveform, the amplitude of the correction growing with 
the decrease of g. Further, the numerically found wave function 4>z(x) is accurately approximated by Eq. ([9| at 
g = 2/3, 1, 2 for A — 0.01, and at g — 2 for A = 0.1, the agreement being qualitative in other cases. 

Comparing these results to analytical solutions (|5| and Q obtained in the uniform system, we conclude that, at 
g = —1 and g = for A = 0.01, at g = — 1, 0, 2/3, 1 for A = 0.1, and all values of g for A = 1, both components in the 
numerically found solutions have approximately equal amplitudes, the difference from the uniform system being that 
the periodic change of the sign of f2(x) gives rise to the array of kinks in <fi-z{x). Another essential deviation from the 
uniform system is observed in the transition between the states with nearly equal and strongly different amplitudes 
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of the two components: as shown above, in the uniform system the transition point, g = 1, does not depend on the 
linear-coupling coefficient, A, while the numerical results demonstrate a strong dependence of the transition on A: in 
the case of the weak coupling, A = 0.01, the transition occurs around g — 0, where the GS of the uniform system 
would remain symmetric; in the case of the moderate coupling, A — 0.1, the transition takes place around g = 2/3, 
i.e., close to the point g = 1 predicted by the analysis of the uniform system; finally, the strong linear coupling, A = 1, 
keeps the amplitudes of the components virtually equal even at g = 2, where the uniform system would feature the 
strong asymmetry, as per solution Q. 



A=0.01 
g=-l 






A=0.01 

g=2/3 





o.4- A=0.1 

g=2/3 



o.o- 



-0.4- 









Figure 1: (Color online) Stationary wave functions (green) and <f>2{x) (red) in the case of the odd modulation function, 

Q.(x) = A sin (ax) (blue), see Eq. pj. The left, middle and right columns correspond to the weak, moderate, and strong linear 
coupling, respectively, with A — 0.01, A — 0.1, and A = 1. The rows running from top to bottom correspond the following 
values of the XPM coefficient: g = -1, 0, 2/3, 1, 2. 
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B. The number of kinks in the ground states 



The number of kinks is the most essential overall (topological) characteristic of the GS [we stress that the kinks 
are always formed in one component, which is denoted (f>^(x), while the other component, (f>i(x), does not change its 
sign]. As said above, in the case of the odd/even spatial modulation of the linear coupling in Eqs. ^ and ([2]), the 
number of kinks is always odd/even too. Analyzing numerical data, we counted those kinks which showed, in their 
core areas, the amplitude of <f>2 that would be no smaller than 5% of the largest value of ip\ . In terms of the present 
analysis, this actually means that the kinks were taken into account if the respective amplitude of 4>2(x) would not 
fall below = 0.02. 

The number of kinks in the GS is presented, as a function of the linear- coupling strength A and modulation 
wavenumber a, in maps collected in Fig. [2j at four fixed values of the XPM coefficient g. The case of g = — 1 is not 
included, as it always gives rise to a single kink, cf. Fig. [T] The latter observation is explained by the fact that the 
attraction between the two components, in the case of g = —1, impedes 4>2(x) to cross zero while 4>\{x) is not too 
small, except for at the central point. On the other hand, at g > the repulsion between the components facilitates 
the appearance of the zero crossings of 4>2(x), thus leading to the increase in the number of kinks. We note that the 
maps pertaining to different positive values of g, i.e., g — 2/3, 1, and 2, are not drastically different. 

Quite naturally, the number of the kinks increases with the modulation wavenumber, a, as the kinks are associated 
with points where Q(x) changes its sign. The same argument explains why the number of kinks increases too with 
the modulation strength, A. 
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Figure 2: The number of kinks in the numerically found ground states, as a function of amplitude A and normalized wavenumber 
a/ao (with ao = 0.09248) of the odd spatial modulation of the linear coupling, see Eq. pi. The maps are displayed at four 
fixed values of the XPM coefficients, g = 0, 2/3, 1, 2. The values of A are shown on the logarithmic scale, to comprise the cases 
of weak, moderate, and strong linear coupling. In this case of the odd modulation, the number of kinks is odd too, see the text. 



Some data concerning the number of kinks has also been collected for the case of the even modulation of the local 
modulation, with function cos (ax) in Eq. pj). The respective map is displayed in Fig. [3] for g = 0. The absence of 
the kink at the central point makes their overall number smaller than in the case of the odd modulation. In this case, 
general features exhibited by the dependence of the number of kinks on parameters A, a and g are qualitatively the 
same as in the case of the odd modulation. 
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1 [ ; 

0.001 0.01 0.1 1 10 

Amplitude 

Figure 3: The same as in Fig. [3] (again, with Qo = 0.09248), but for the case of the even modulation function, cos (cm;), in Eq. 
Q, and g — 0. The number of kinks is even in the case of the even modulation, see the text. 

IV. CONCLUSION 

The objective of this work was to identify the GSs (ground states) in two-component models based on the system 
of two NLSEs/GPEs coupled by the linear interconversion and nonlinear (XPM) interaction, with the self-repulsion 
acting in each component. A well-known property of such systems is that their GSs may be represented by either 
symmetric or asymmetric states, depending on the strengths of linear and XPM couplings. In the present model, 
we have extended the idea of Ref. [21] and considered in detail the periodic odd, or even spatial modulation of the 
linear coupling, f2(x). The repeated change of the sign by O(x) produces a strong impact on the structure of the 
GS in the two-component system, as the relative sign of the components tends to be locked to the local sign of fl. 
Accordingly, the GS of the trapped system features one or several kinks in one component, while the other component 
does not change its sign. The parity of the number of kinks coincides with that of the modulation function. The 
interplay of the modulation of fl(x) with the trapping potential also strongly affects the transition between the states 
with nearly equal and strongly different amplitudes of the two components. These results were obtained by means 
of the systematic numerical analysis, and also in the analytic perturbative form. The model may be realized in 
nonlinear optics, in terms of periodically modulated slab waveguides. Also, it may also be implemented in BECs, 
when two-photon Raman transitions are used. 

A challenging problem is to extend the analysis to two-dimensional binary systems, with a two-dimensional 
(checkerboard-patterned) modulation of the linear coupling. This two-dimensional setting may also be realized in 
the photonic-crystal medium. In 2D, as as already discussed in [21], it should be possible to generate vortices and 
(in more complex settings) other types of topological defects. In fact, this method of vortex generation is completely 
analogue to the methods of creating "artificial" magnetic fields employing spatially dependent Berry's phase |31| . 
or Raman couplings [32] • Another interesting extension of the present studies is to perform a similar analysis for 
quantum systems, formulated in terms of an appropriate Bose-Hubbard model, i.e. in the strongly correlated limit, 
where one expects the appearance of quantum solitons (cf. |33j and references therein) . 
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